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Abstract 

We study a simplified version of tlie Standard Electroweak Model and intro- 
duce the concept of the physical gauge invariant effective potential in terms of 
matrix elements of the Hamiltonian in physical states. This procedure allows 
an unambiguous identification of the symmetry breaking order parameter and 
the resulting effective potential as the energy in a constrained state. We ex- 
plicitly compute the physical effective potential at one loop order and improve 
it using the RG. This construction allows us to extract a reliable, gauge in- 
variant bound on the Higgs mass by unambiguously obtaining the scale at 
which new physics should emerge to preclude vacuum instability. Compari- 
son is made with popular gauge fixing procedures and an "error" estimate is 
provided between the Landau gauge fixed and the gauge invariant results. 
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I. INTRODUCTION AND MOTIVATION 



Considerable effort has been invested in understanding the Higgs sector of the Standard 
Model. Reliable constraints on the Higgs boson mass are important in determining the 
energy scales for collider experiments that probe the electroweak sjTiimetry breaking sector. 
Alternately, should the Higgs be discovered, its properties will help elucidate high-scale 
physics, (see ^ for an early review). In the Standard Model the requirement that the 
conventional effective potential Kfr(v^) have its global minimum at the electroweak scale has 
been used to obtain a relation between a lower bound on the mass of the Higgs boson and 
the scale at which new physics should appear (for recent reviews see 0,0). That scale is 
related to the value of (p (expectation value of the Higgs field) at which Ves develops a new 
deeper minimum (which depends on the mass of the Higgs). Recently Loinaz and Willey 
have pointed out a difficulty with this procedure. When the contributions from the gauge 
sector of the electroweak theory are included in the effective potential, the value of which 
minimizes Ves is gauge dependent. The gauge dependence of the effective potential was 
already recognized by Dolan and Jackiw in their early formulation of effective potentials . 
Although in some specific gauges the contribution of the gauge sector may be a perturbative 
correction to the contributions from the scalar plus heavy (top) fermion- Yukawa sector, the 
gauge dependence implies that no error estimate can be made and raises questions on the 
reliability of such bounds. The conventional V^s, defined in terms of IPI Green's functions at 
zero momentum is an off-shell quantity and in a gauge-fixed formulation is inherently gauge 
variant and therefore not uniquely defined. Using the effective potential to study stability 
or metastability implicitly assumes that Ves is associated with the energy of (space-time 
constant) field configurations. In scalar field theories the effective potential is proven to be 
the energy of a constrained state , but such proof is lacking in gauge theories. Since the 
effective potential as calculated in gauge- fixed formulations is explicitly gauge dependent, it 
cannot generally be identified with the expectation value of the physical Hamiltonian in a 
physical state. There are, however, gauge independent quantities that can be extracted from 
the effective potential. Nielsen identities ||T^ have been used to prove that the difference of 
the values of the (gauge fixed) effective potential evaluated at extrema are gauge invariant, 
as well as the nucleation rate for bubbles in a first order transition when calculated from 
the effective potential |TT|. However providing a lower bound on the Higgs mass requires 
to estimate values of the expectation value of the scalar order parameter, which is a gauge 
variant quantity even at the minima of the effective potential. 

Thus it is important to provide a formulation in which the effective potential is a gauge 
invariant function of a gauge invariant order parameter and can be interpreted as an en- 
ergy function. There have been efforts to formulate a gauge invariant effective action |T^ 



and consequently an effective potential, but the formalism involved is formidable and its 
calculational implementability rather unwieldy. 

Recently a formulation that allows to obtain a gauge invariant effective potential as the 
expectation value of the Hamiltonian in physical states has been developed within a different 
framework |T^. We refer to this as the Physical Effective Potential (PEP). In this article 



we apply the formulation proposed in to solve the problem of the gauge dependence 
in a slight variant of the model studied in Although this model is a simplified Abelian 
version of the Standard Model, it serves to demonstrate the utility of the gauge-invariant 
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effective potential in providing a gauge-invariant estimate of a vacuum instability scale. 

This article is organized as follows: in section II we present the model to be studied 
and the relevant aspects of the gauge invariant formulation provided in [|13| and adapted 
to include the fermionic sector. In this section we define the physical (gauge invariant) 
observables including the order parameter that provides an unambiguous signal for symme- 
try breaking. In section III we explicitly construct the one loop effective potential as the 
expectation value of the physical Hamiltonian in gauge invariant states constrained to give 
a space-time constant expectation value of the gauge invariant order parameter. We also 
provide the MS renormalization of this effective potential. In section IV we compare our 
results to those obtained from the gauge fixed formulation in general R^^u gauges. Section V 
is devoted to a RG improvement of the gauge invariant effective potential and to an unam- 
biguous determination of the lower bound on the Higgs mass in this model, providing "error 
estimates" for the quantities obtained in the gauge variant formulation. In Section VI we 
present some brief numerical results. 

Section VII summarizes our conclusions and suggests possibles avenues to extend the 
gauge invariant construction to non-Abelian gauge theories. 

II. THE GAUGE INVARIANT DESCRIPTION: 

The focus of our study is the Abelian Higgs model with an axial coupling of the gauge 
fields to fermions. The Lagrangian density is 



C = -^F^^'F^, + Dy^D^<f) - /iVV - A(0V)' + ^ [^P+ y (01 + ^7'02) ^] (2.1) 
D,<P = id,<P + ieA,<P) ; p=r{d,~ ^el'^) ; = i= (0i + z^^) (2.2) 



Our goal is to define and explicitly evaluate the gauge invariant effective potential, in 
terms of a gauge invariant order parameter. The definition of the effective potential that 
we use is in terms of the expectation value of the physical Hamiltonian in physical states 
constrained to provide a fixed expectation value of the physical order parameter. 

The program for the construction of such an effective potential requires i) the identifica- 
tion of the physical (gauge invariant) states of the theory, ii) the identification of an order 
parameter that is invariant under the local gauge transformation but transforms nontrivially 
under the rigid global symmetry transformation and thus gives information on spontaneous 
symmetry breaking iii) construction of the gauge invariant effective potential as the expec- 
tation value of the Hamiltonian in constrained physical states. Since the concepts behind 
the construction are not part of the standard lore, we highlight below the most relevant 



aspects of the formulation, for more details see |T3 . 

Such a description is best achieved within the canonical formulation, which begins with 
the identification of canonical field variables and constraints. These will determine the 
classical physical phase space and, at the quantum level, the physical Hilbert space. 

The canonical momenta conjugate to the scalar and vector fields are given by 

n° = ; W = A' + VM° = -E' (2.3) 

Tft = + ze^V ; 77 = 0^ - ieA^(f)^ (2.4) 
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The Hamiltonian is therefore 



d^a; <{ ■ n + ttV + (V0 - ieA^) ■ (V0^ + ieAcp^) + ^(V x Af+ 



H(j)^(f) + A(0V)^ + ^0 V ■ n - ie{TT(j) - vrV^) + i^^^^l^^ 



+ 



(2.5) 



where a ; /5 are the Dirac matrices. 

There are two main methods for quantizing gauge theories, the first one originally due to 
Dirac [Q(see also |T^) begins by identifying the first class (mutually commuting) constraints 
and projects the physical states by requiring that these are simultaneously annihilated by all 
first class constraints. Physical operators commute with all of the first class constraints. The 
second, and most used method, "fixes a gauge" , converting the set of first class constraints 
into second class constraints (non-commuting) and introducing the Dirac brackets. This is 
the popular method of dealing with the constraints and leads to the usual gauge-fixed path 
integral representation [0] in terms of Faddeev- Popov determinants and ghosts. Although 
this second method is the most popular as it is easily translated into a path integral language, 
it has the drawback that the physical quantities are more difficult to extract, and though 
S-matrix elements are gauge invariant, off-shell quantities generally are not. 

In order to avoid ambiguities and to define a physical order parameter and effective 
potential (an off-shell quantity) we choose to use the first method. 



In Dirac's method of quantization |14,15| there are two first class constraints which are: 



and Gauss' law: 



g{x,t) = v-n 
p 



6C 

Ja^ 



p = 



ie{n(f) — 7r^0^) — i-ip'^'j^ip 



(2.6) 

(2.7) 
(2.8) 



with p being the matter (complex scalar and fermionic) field charge density. 

Gauss' law can be seen to be a constraint in two ways: either because it cannot be 
obtained as a Hamiltonian equation of motion, or because in Dirac's formalism, it is the 
secondary (first class) constraint obtained by requiring that the primary constraint ( p.6|) 
remain constant in time. Quantization is now achieved by imposing the canonical equal- 
time commutation relations 



U%x,t),A%y,t) =-i5{x-y) ; W{x,t), A^ {y,t) = -r5^^{x - y) 



(2.9) 



along with the usual canonical commutators for the scalar field and its canonical momentum 
and anticommutators for the fermionic fields. 

In Dirac's formulation, the projection onto the gauge invariant subspace of the full Hilbert 
space is achieved by imposing the first class constraints onto the states. Physical operators 
are those that commute with the first class constraints. Since 11° (x, t) and Q{x, t) are gener- 
ators of local gauge transformations, operators that commute with the first class constraints 
are gauge invariant ||13 
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Following the steps of reference |T3[ we find that the fields and canonical conjugate 
momenta 



U(x) 



(x) exp 



TT[x) exp 



ie J d^yA{y) -VyGiy-x) 
-ie I d^yA{y)-VyG{y-x) 



1 

71 



[^1 + 1^2) 



'^{x) = exp 



d'yA{y)-VyG{y-x) 



ip{x) 



with G{y — x) the Coulomb Green's function that satisfies 

V'^G{y-x) =6%y-x) 



(2.10) 
(2.11) 
(2.12) 

(2.13) 



are invariant under the gauge transformations |T3| . Furthermore writing the gauge field into 
transverse and longitudinal components as follows 



A{x) = Al{x) +At{x) 

V X Al(x) = ; V ■ At(x) = 



(2.14) 
(2.15) 



it is clear that the "transverse component" At{x) is also a gauge invariant operator. The 
fields At ; $ ; \& and their canonical momenta are gauge invariant as they commute with the 
constraints. The momentum canonical to A ; 11 is written in terms of "longitudinal" and 
"transverse" components 



U(x) =Ul(x) + Ut(x) 



(2.16) 



Both components are gauge invariant. 

In the physical subspace of gauge invariant wave-functionals, matrix elements of V ■ 11 
can be replaced by matrix elements of the charge density p, since matrix elements of Gauss' 
law between these states vanish. Therefore in all matrix elements between gauge invariant 
states (or functionals) one can replace 



Ilr(x) 



P 



(2.17) 



with the charge density (a gauge invariant operator) written in terms of the gauge invariant 
fields as 



P 



%e 



($(^)n(^) - ^\m\y)) - ^^(y)7'^ 



(2.18) 



This procedure is tantamount to solving the constraints in the physical space |T5 
Finally in the gauge invariant subspace the Hamiltonian becomes 

11 = ^ S-x I^IIt ■ fir + n^n + (V$ - itAr^) ■ (V$"f + ieAr$^) + ^(V x At) 



-la 



- d-^y d^xp{x)G{x - y)p{y) 



(2.19) 
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Clearly the Hamiltonian is gauge invariant, and it manifestly has the global U{1) (chiral) 
symmetry under which 



$(f) ^ $(f)e^'^ ; ^(f) e~'^ ^ ^(f) (2.20) 

with ip a constant real phase, 11 transforms with the opposite phase and Aj- is invariant. 

The Hamiltonian written in terms of gauge invariant field operators ( |2.19| ) is reminiscent 
of the Coulomb gauge Hamiltonian, but we emphasize that we have not imposed any gauge 
fixing condition. The formulation is fully gauge invariant, written in terms of operators that 
commute with the generators of gauge transformations and states that are invariant under 
these transformations. The similarity to the Coulomb gauge Hamiltonian is a consequence of 
the fact that in this Abelian theory. Coulomb gauge displays explicitly the physical degrees 
of freedom. 

There is a definite advantage in this gauge invariant formulation: the (composite) field 
$(x) is a candidate for a locally gauge invariant order parameter. The point to stress is the 
following: this operator is invariant under local gauge transformations but it transforms as 
a charged operator under the global gauge transformations generated by Q = / d^xp{x), that 
is 

e^°'?$(f)e-*"'3 = e''"<!>{x). (2.21) 
Because the gauge constraints annihilate the physical states and these constraints are the 



generators of local gauge transformations |13,0|, these states are invariant under the local 



gauge transformations and any operator that is not invariant under these local transforma- 
tions must have zero expectation value. The local gauge symmetry cannot be spontaneously 
broken; this result is widely known in lattice gauge theory as Elitzur's theorem [T^. How- 



ever, the global symmetry generated by the charge Q can be spontaneously broken and the 
expectation value of a charged field signals this breakdown. 

From this discussion we clearly see that a trustworthy order parameter must be invariant 
under the local gauge transformations, thus commuting with the gauge constraints, but must 
transform non-trivially under the global gauge transformation generated by the charge. The 
field $ fulfills these criteria and is the natural candidate for an order parameter. 

At this stage, having recognized the physical states one could prefer to pass to a path 
integral representation of the vacuum-in to vacuum-out transition amplitude. This can 
now be done unambiguously by carrying out the usual procedure in terms of phase space 
path integrals with the gauge invariant measure VUt'DAt'D^V^'^VIIVW ■ ■ ■. There is no 
need for "gauge fixing". In the resulting action (of the form pq — H), the instantaneous 
Coulomb interaction can be re-written by introducing an auxiliary field, and the integral 
over the canonical momenta can be carried out explicitly leading to a Lagrangian form. The 
resulting Lagrangian leads to Feynman rules that are very similar to those in Coulomb gauge 
and allow the perturbative calculation of wave function renormalization constants needed 
below. The gauge invariant effective potential can also be computed in this path integral 
representation, but we prefer to provide its explicit construction from the Hamiltonian as 
such construction displays more clearly the identification of the effective potential as the 
energy of a constrained state. 
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III. THE EFFECTIVE POTENTIAL 



We are now in position to define the gauge invariant effective potential. Consider the 
class of gauge invariant states |T,x) characterized by the condition that the expectation 
value of the gauge invariant order parameter in this state is nonzero and space-time 
constant 

{T,xmx)\T,x) 1 . ^. .... 

:{Xi + iX2) = X (3.1) 



(T,x|T,x) V2' 



In this notation T indexes the states within the set characterized by eq. |3.1| . The effective 
potential is defined as the minimum of the expectation value of the Hamiltonian density 
over this class of constrained states [pHH], namely 



with H being the gauge invariant Hamiltonian given by equation (|2.19|) and Q the spatial 
volume 10. By construction in terms of gauge invariant states (or functionals) and the 
gauge invariant Hamiltonian, this effective potential is gauge invariant. The minima of this 
effective potential are obtained from Kfr(x) by further minimizing with respect to Xi.2- 
It is convenient to separate the expectation value of $ as 

'^(^) = ^(^1 + *^2) + ^(^i(^) + *^2(x)) (3.3) 
H(x) = i=(Hi(x)-zH2(x)) (3.4) 

The one-loop correction (formally of 0{h)) to the effective potential is obtained by keeping 
the linear and quadratic terms in rii^2 in the Hamiltonian, however the linear terms will not 
contribute to the effective potential because their contribution vanishes upon taking the 
expectation value in the state \T,x)- Thus keeping only the quadratic terms in 7712 we 
obtain 

H, = n(C\x\' + ^\xA+H^ + H' 



H^ = \j d'x {n| + (V X Af + e'\x?Al + H? + H^ + (Vr/i)^ + (Vr^2)'+ 

r^?(;.2^3A|xP)+%'(/i' + A|xP)}-^lxP/ j d^yIi,{x)G{x - m2{y) 

= j d^x[^^[~ia-W + y\x\P]-^} (3.5) 

where we have performed a rigid chiral phase rotation under which the Hamiltonian is 
invariant. The transverse components At describe a field with mass m\ = e^|xP and only 
two polarizations. The fermionic part is recognized as a free Dirac fermion with mass Mp = 
y\x\- The fermionic Hamiltonian can be diagonalized in terms of particle and antiparticle 
creation and annihilation operators of the usual form, and upon using their anticommutation 
relation we obtain 
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H^ = T.i: {^^(k) [blA,^ + 4, } - 2 E ^^(k) 

k a=l,2 k 

u^{k) = v/F + M|[|x|] ; Mp[x] = y\x\ 



(3.6) 
(3.7) 



Following we write the bosonic Hamiltonian (H^) in terms of the spatial Fourier 



transform of the fields and their canonical momenta in terms of which the quadratic part of 
the Hamiltonian finally becomes 



U,ik)U,i-k) + cuUk)mik)vii-k) + U2ik)U2i-k) 



uj^{k) 
A;2 



+ V2ik)v2{-k)co'Jk)\ (3.8) 



where the frequencies are given in terms of the effective masses as 



,^^{k) = e + M^[x] 

cuUk) = k' + MUx] 

u^ik) = e + M'^lx] 



UJ 
.2 



2| |2 



M^lx] 
M]j[x]=fi' + 3X\x\ 
M'[x]=fi' + X\x\'. 



(3.9) 
(3.10) 
(3.11) 



The last two terms can be brought to a canonical form by a Bogoliubov transformation. 
Define the new canonical coordinate Q and conjugate momentum P as 



Il2{k) 



k 



ujrik) 



P{k) 



V2ik) ^ ^Qik) 



(3.12) 



in terms of which the last term of the Hamiltonian (|3.8| ) becomes a canonical quadratic form 
with the plasma frequency 



ulik) = uj^g{k) 



uj^{k) 



[e + M'Jx]][k' + MM]/k' 



(3.13) 



There are four physical degrees of freedom. The modes with frequency uxik) are the two 
transverse degrees of freedom, the mode with frequency Unik) is identified with the Higgs 
mode. In absence of electromagnetic interactions (e = 0) the mode with frequency ujp{k) 
represents the Goldstone mode whereas in equilibrium, namely at the minimum of the tree 
level potential, when 



Ixl = y — /i^/A, it represents the plasma mode which is identified as 
the screened Coulomb interaction, and the transverse and plasma modes all share the same 
mass. A detailed discussion of the dispersion relation of the bosonic excitations has been 
provided in [ll3| . 

The quadratic Hamiltonian is now diagonalized in terms of creation and destruction 
operators for the quanta of each harmonic oscillator. The ground state is the vacuum 
for each oscillator and is the state of lowest energy compatible with the constraint ( ^.1[ ). 
Therefore the one loop {0{h)) contribution to the effective potential is obtained from the 
zero point energy of the bosonic oscillators and the (negative) contribution from the "Dirac 
sea" given by the last term in (|3.6| ). Therefore accounting for the two polarizations of the 
transverse components we find: 
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I r d^h 

Kff(lxl) = + - y (^[2^t(A:) + i^H{k) + i.,{k) - Au;^{k)] (3.14) 

The normalized state that satisfies (|3.1| ) and gives the minimum expectation value of the 
Hamiltonian, thus determining effective potential via ( p.2| ) is given by 



Ix) = \x)t, ® \x)t, ® \x)h ® \x)p ® \x)i 



(3.15) 



i.e. a tensor product of the harmonic oscillator ground states for the two polarizations ( 
Ti 2 ), the Higgs mode (H), the "plasma" mode (p) and the fermionic Fock state (F) (the 
"Dirac sea"). The functional representation for the bosonic Fock states is simply a Gaussian 
wave-functional. This state hy construction is gauge invariant. 

The k-integrals in the final form of the gauge invariant effective potential ( |3.14D 
are carried out in dimensional regularization by the replacements: Jd^k/{27i)^ — »• 
/ d'^~^k/{2TrY~^; = d, and 2 — > — 2 transverse modes, d = 4 — e (see the appendix) 
and we obtain the following (unrenormalized expression) one loop effective potential 



Kff(x) = V4x) + 



1 



(47r)2 



-A 
( In 



r 



Mhlx] 



In 



a 

Ml[x] 



+ 



+ 



{MT[x]+MMf 



3 



--Mt[x]M,[x](Mt[x] - M,[x]f - MUx] In 



Mlix] 



3 



with 



A. 



A-d 



7 + ln[47r] 



(3.16) 



(3.17) 



is the standard MS UV subtraction in d-dimensions, k is the renormalization scale and 7 is 
the Euler-Mascheroni constant. 



A. MS Renormalization 



The theory is renormalized by coupling, mass and wavefunction renormalizations. 

66 = Z^tR ; A;, = Zx\r ; Ub = Zy-yn ; i4 = ^^2/^^ (3.18) 



= ZlAn ; $fe = Zl*^ ; = zI^r (3.19) 

The gauge invariant quantization ensures that all the Ward identities associated with 
the gauge symmetry are now trivially satisfied. Here we work in the MS renormalization 
scheme and write for the various couplings and wavefunction renormalizations Z = 1 + 6Z, 
expanding 6Z in powers of the couplings and absorbing only the MS divergences proportional 
to Ag. There are two equivalent procedures to renormalize the effective potential: 
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1. The terms proportional to in ( p.l6| ) are absorbed in a partial renormalization of 
couplings and /z^. This renormalization, however, does not render finite the scattering 
amplitudes. The latter require a further renormalization by the proper wave function 
renormalization constants, which at the level of the effective potential is absorbed 
in a renormalization of xl ~^ ^^Xji- The wave-function renormalizations must be 
calculated differently from the effective potential, since their calculation requires one- 
loop contributions at non-zero momentum. 

2. After restricting the Hamiltonian to the gauge invariant subspace, one can pass on to 
the Lagrangian density in terms of gauge invariant variables, and write it in terms of 
the renormalized couplings, masses and fields plus the counterterm Lagrangian. The 
counterterms are then required to cancel the divergences in the proper IPI Green's 
functions. It is at this stage that we use the path integral representation discussed 
previously to compute the wave function renormalization constants from the one loop 
self-energies and vertex corrections. Such a calculation, although straightforward, 
involves non-covariant loop integrals which are performed with the help of the appendix 
in reference 



T9|. 



For the renormalization of the effective potential only scalar wave function renormalization 
is needed beyond the cancellation of the terms proportional to in (|3.16| ). A lengthy but 
straightforward calculation of the relevant one loop diagrams provides the renormalization 
constants listed in Appendix B. 

In terms of the renormalized couplings and expectation value of the gauge invariant 
scalar field, the gauge invariant effective potential in the MS renormalization scheme is 
given by (here we drop the subscripts (R) in the renormalized quantities to avoid cluttering 
of notation, but all quantities below are renormalized) 



(4vr)2 [ 2 \ 
i(M|[x]-M|[x])'fln 



In 



MM 



+ 



{MT[x] + MMf 



-Mt[x]MM{Mt[x] - Mg[x]Y - M^M In 




MUx] 



+ 



(3.20) 



IV. COMPARISON WITH GAUGE-FIXED RESULTS: 

A comparison with the effective potential in covariant R^^^u gauges is established by 
adding a gauge fixing and Faddeev- Popov terms to the Lagrangian density eq.( pJ]) : 

C + Cgf + Cfp (4.1) 

C,f = -^^{d,A^ + iuecl^^f (4.2) 

Cpp = —cd'^c — ^ue'^cc(j)i (4.3) 

Special cases of this covariant gauge fixing procedure include the Landau gauges (^ = 0), 



't Hooft gauges {u = (0i)q, the tree-level Higgs vev) and Fermi gauges {u = 0) 
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where we have chosen the symmetry breaking expectation value along the 0i direction for 
convenience. At this stage we can just quote the results of reference adapted to the 
case treated here of an axial vector coupling with Nf = 1, for arbitrary number of flavors, 
the one-loop contribution from the fermion will be multiplied by A'^^. With (f = the 
expectation value obtained in the gauge-fixed path integral, we obtain the one loop effective 
potential in the MS scheme after renormalization of /x^, couplings and wavefunctions |Q: 



r . 1 1 2 2 A 4 IH^h] 



3B^ 
+- 



+ 



4 (47r) 



In 



4 (Airf 



In- 



(4vr) 



In 



y2 



In 



2G^ 



4 (47r) 
Ik* \^ 



3 
2 

In 
In- 



4 {Anf 



3 

" 2 

-A, 

(4vr) 



(4.4) 



where 



kl 



I? + 3A(^^ 

2 2 

e if 



± 



-v/(^2 + ^^2) (^2 + ^^2 + 4^e2^ _ ^)) 
2 2 



(4.5) 
(4.6) 
(4.7) 



(4.8) 
(4.9) 



H, B, Y and G denote contributions from Higgs, vector boson, fermion and Fadeev-Popov 
ghost loops, respectively. 

Although the "effective masses" H[(f] ; B[(f] ; Y[(f] are formally similar to the Higgs, 
transverse and fermion effective masses given in equations ( ^ , ^ , 3.101) upon the replace- 
ment (y9 — >- Ixl, we want to emphasize that ipex is the expectation value of a gauge dependent 
field in a gauge-fixed state, whereas x is a true gauge invariant order parameter, the ex- 
pectation value of a gauge invariant field in a gauge invariant state. The UV pole in the 
gauge-fixed effective potential ([4.4|) which is not removed by couplings and wave-function 
renormalizations has been discussed in detail in . It can be removed by a shift in (p consis- 
tently in the loop expansion and is a consequence of the particular choice of R^^u gauge- fixing 
that breaks the global U{1) symmetry explicitly. 

Even for the Landau gauge, corresponding to the choice ^ = 0, leading to k^ 
Xip"^; k^ = the difference of the gauge sector contributions to the gauge fixed and the 
gauge invariant effective potential is clear. 

The first physical quantity that we must compare is the value of both the gauge invariant 
and the gauge fixed effective potentials at their respective minima. The minima are obtained 
from the conditions 
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' ^Kff [1x1 
d\x\ 



dip 







Writing the stationary values of |x| ; 93 in a formal loop expansion as 

Ixlext = Ixlo + Ixll H ; V'ext = v^o + H 

we find the following relations that fix the minima to one-loop order 



dK'ff [\X\ 



\X\i 



d\x\ 

^eff 



Ixlo 



mf [\x\ 



d\x\ 



Ixlo 



d'K'ff [\x\ 



d\x\ 



Ixlo. 



dip 
dip 



^0 



-1 



(4.10) 

(4.11) 

(4.12) 
(4.13) 
(4.14) 

(4.15) 



fa , 



the solutions of eq.( [4.12| , [4.14|) are obviously the tree level vacuum expectation values |x|o = 
— /i^/ A but the solution of eq. (^4.15|) is dependent on the gauge fixing parameters ^, u 



^0 = 

1^. However, upon inserting the solutions of the minima equation consistently up to one 
loop in the respective expressions for the effective potentials, we find that the value of the 
minima for the gauge-fixed and the gauge invariant effective potential are the same. This 
is an important result: whereas the gauge invariant effective potential as constructed has 
the meaning of an energy of constrained physical states, no such interpretation is available 
for the gauge-fixed result. However at the minimum, we see that the gauge- fixed effective 
potential agrees with the gauge invariant result and therefore the values of the extrema 
provide gauge invariant information on the energetics of constrained states. 

The second derivative of the gauge fixed effective potential at the minimum is seen to 
be a gauge dependent quantity. Although the pole mass of the excitations are on-shell 
quantities and must therefore be gauge invariant, the second derivative of the effective po- 
tential corresponds to the value of the one-particle Green's function at zero four momentum 
transfer. 



dip^ 



2Xipl + 



J ip=ipo 



16772 



6X' 



4 In 



2Xipl 



+ 9e^ In 



+ 



2X^ + lx^e'+{^e'y 



\J-^-X^e'-16N,yt 



31n^-l 



(4.16) 



We note that this expression is divergent in the Landau gauge (^ = 0). Direct calculation 
in Landau gauge shows that the divergence is due to the Goldstone boson term in the effective 
potential. 
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V. RG IMPROVEMENT AND THE HIGGS MASS BOUND 



The vacuum stability bound on the Higgs mass arises from imposing the requirement 
that the electroweak vacuum be the global minimum of the effective potential. In fact, the 
Standard Model is thought to be the low-energy effective theory of some more fundamental 
high-scale theory. Thus it is only consistent to demand that the electroweak vacuum is 
the absolute minimum up to the scale at which the effects of the 'new physics' (those not 
incorporated into the low-energy effective theory) become significant. In the context of the 
effective potential we associate this scale with a value of the elementary scalar field and 
hence insist that the electroweak minimum is the global minimum of the effective potential 
up to that value of the field (however, see [0). For the elementary scalar field effective 
potential, however, the value that the effective potential assumes at some expectation value 
of the scalar field is explicitly gauge dependent. Thus, the statement that the value of 
the effective potential exceeds that of the electroweak minimum up to some scale A is also 
gauge dependent. This gauge dependence in turn infects the lower bound on the Higgs mass 
derived from that ansatze. 

The PEP does not suffer from this gauge ambiguity. The condition that the electroweak 
minimum is the global minimum up to some maximum value of the gauge invariant order 
parameter provides a gauge invariant means of defining the vacuum instability scale and 
through that a gauge invariant lower bound on the Higgs mass. 

Using the PEP, we demonstrate how, to derive the lower bound on the Higgs mass in the 
toy model, which displays the same qualitative features as the Standard Model with respect 
to vacuum stability. Unlike the conventional approach using the elementary field effective 
potential, the bound obtained from the PEP formalism will be manifestly gauge invariant. 
The extension of the approach to nonabelian gauge theories such as the Standard Model 
requires the extension of the PEP formalism to that context, we see no serious obstacles to 
such a program. 

It is known (in the context of the elementary field effective potential) that the simple 
perturbative effective potential is often inadequate for the study of the theory at large 
field values due to large logarithms which ruin the convergence of the perturbative loop 
expansion |lT7|,p|. The range of validity of the perturbative expansion can be enhanced, 



however, through the use of the renormalization group. The invariance of the full (all- 
orders) effective potential under a change in MS renormalization scale can be written as a 
first-order differential equation 

K^Kff (V5, yu, A, e, yt, ^, k) = (5.1) 

where ^ denotes a generic gauge parameter, which always arises in defining the gauge-fixed 
elementary field effective potential. This equation may be solved using familiar methods in 
terms of the /3 and 7 functions. These in turn may be calculated in a loop expansion from 
the MS counterterms of the theory. Whereas the unimproved effective potential was reliable 
only for field values at which 

(6.2) 



IGtt^ k 
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the solution to the RG equation will be reliable as long as the running couplings (generically 
g'^{s)) are small. 

Similar considerations are applicable to the PEP (but now without the gauge parameter 
^). The PEP is also independent of changes in renormalization scale. The PEP has been 
identified as a matrix element of the physical Hamiltonian in (physical) states. The Hamil- 
tonian is part of the energy momentum tensor which is a conserved current. Conserved 
currents do not acquire anomalous dimensions and are finite after field independent sub- 
tractions (normal ordering in the free field vacuum) in terms of the renormalized parameters 
and fields. In dimensional regularization these normal ordering divergences vanish identi- 
cally. Alternatively, that the PEP is renormalization group invariant can also be seen in the 
same manner as for the elementary field effective potential, by returning to the path inte- 
gral, expressing the PEP as a sum of $"rSr^ at zero external momentum via the momentum 
expansion of the corresponding effective action. The multiplicative renormalization factors 
cancel, and since the renormalized effective potential equal to the bare effective potential it 
is also K -independent. The analogous RG equation may then be written: 



«:^Kfr (x,/^,fi', «:) 



d ^ d d d 



VeS iX, Ai, g,fi) =0 



(5.3) 



where 



Pa 



K- 



dg_ 
K d^ 



1 K dZ^ 
^ dn 2 Z<j, dn 
Kdfi Ik dZ^2 
^dn 2 Z^2 dn 



7/. 



(5.4) 



or, writing $ = s$i (where $.j is some arbitrary initial scale, for example the electroweak 
scale) 



Ok ^ do 



d_ 



d_ 

ds 



(5.5) 



Here g represents the set of couplings X, g, and yt- Using dimensional analysis we can also 
write 



d d d ^ 

OK OjJ, OS 



0. 



Combining eq.B3| and eq.|^ gives the RG equation 



(5.6) 



~'~ dg 



(7* + l)^^ + rf 



Vcs {sxi,fi,g,n) = 0. 



(5.7) 



It can then be shown that solution of the MS RG equation for the effective potential is 

VeS (sXi.gi, k) = [Cis)f VeS (Xi, g ("S, ^i) , /i (s, /ij) , k) (5.8) 
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where 



C(s) = exp 

g{s,gi) is the solution to the equation 

dg{s) 
ds 

m = 9 



log s 



7# (x) + 1 



-dx 



and 



log s 



7$(x) + 1 



dx 



(5.9) 



(5.10) 



(5.11) 



It is convenient to separate the tree-level Ves and its one-loop corrections by writing eq.|578 
in the form 



i {m\s) + Am\xi)) X- + \Ks{s, X^)x\ 



(5.12) 



Aeff(s,Xi) = A (s) + AA (x») 



(5.13) 



where Am^ and A A contain the loop corrections to 

The /3 and 7 functions and \4ff can then be calculated to the desired loop order. It has 
been shown that using the n-loop effective potential and the n + 1 loop /5 functions sums 
up all logs of the form log s, log^ s, . . . , log" s [0,^. The one-loop MS /5 functions for our 
model are given explicitly in Appendix B. A A can be extracted directly from eq.( p.l6|) and 
to one-loop order is 



AA(s,Xi 



In 




M^AX^ 



In 



^(M|[x.]-M|[x.])'fln 



{MT[X^]+Mg[x^]f 



--MT[x^]MM{MT[x^] ' M,[x^]f - M^Xi 



In 



MUX^ 




(5.14) 



In the Standard Model the top quark term dominates (3x. In this toy model the fermion 
plays the same role. It is the large negative fermion loop contribution that drives A(s) 
negative at large s. At the leading-log level, if the tree-level RG-improved effective potential 
is run with the one-loop (3 function, at large fields it is well-approximated by neglecting 
the term quadratic in x and is given by V^ff[s,Xi] = i^{s)[XiC{.^)]'^- Very near the value 
of s at which A(s) turns negative, Ves falls below the electroweak minimum. We define 
the scale s at which Ves{s) = Kfr(sEw) to be the vacuum instability scale, syi. If the 
electroweak minimum is to remain the global minimum of the effective potential, the theory 
must somehow be incomplete and new high-scale physics must contribute in some important 
way before the vacuum instability scale. 
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To extract a lower bound on the Higgs mass, we choose a vacuum instabihty scale syi 
where by definition Kfr(svi) = Kfr(sEw)- In the approximation that s 1, the contribution 
to Ves{syi) from the terms is small, and Vesisyi) ~ |Aefr(svi)xf [Cl-^vi)]"^- Further, for 
s ^ 1, the value of s at which Ves{s) = is very close to that at which it equals V^fr(sEw)- 



Thus it is numerically a good approximation to take [p^-p6 

Vosisyi, x^) ~ \Ks{sYi)xt [C{sYi)f ~ 0. (5.15) 

Since ({s) is positive, eq. |5.15| implies 

Aeff(svi) = A(svi) + AA(svi) = ^ A(svi) = -AA(svi) (5.16) 

Thus the constraint that the electroweak minimum be the absolute minimum of the effective 
potential up to some high scale translates to a high-scale boundary condition on A(s). Run 
down to the low scale, the result is a value of Aj below which the electroweak vacuum is 
unstable. This can in turn be converted to a lower bound on the Higgs pole mass. 

This is completely analogous to the procedure with the elementary field effective po- 
tential. In that case, however, the expression for AA is explicitly dependent on the gauge 
parameters, which results in a gauge dependence in the Higgs mass bound |Q . In the current 
formulation, AA is by construction gauge- invariant. The steps all follow through in exactly 
the same fashion, but the Higgs mass bound obtained in the end is gauge invariant. 



VI. NUMERICAL RESULTS 

Here we demonstrate the formalism developed in the previous section by applying it to 
obtain numerical values for Aj, assuming values for the other couplings at the electroweak 
scale and assuming some value for syi. Since this is only an abelian toy model, we cannot 
draw conclusions about the Standard Model. Instead, we compare the results of the gauge 
invariant formulation with the elementary field effective potential formulation within this 
model for several different sets of parameters. 

For illustrative purposes it will suffice to use only one-loop P functions and the one- 
loop corrections to the effective potential, which will sum the leading logs. Furthermore, 
at the leading log level it is consistent to run the couplings in the tree-level part of the 
effective potential only and to leave the couplings in the one-loop terms fixed at their initial 
values. The additional effort required to reduce the k dependence of the results and improve 
numerical accuracy by calculating higher loop effects is not justified for a toy model. 

In Figure 1 we plot the log of the vacuum instability scale syi as a function of Aj and 

for = 0.5 and k = f , as calculated using the PEP. We have also performed a similar 



calculation using the Landau gauge elementary field effective potential (setting ^ = in |4.4|) 
and found the results to differ only by a few percent. This is shown explicitly in Figure 
2. Here the difference between the logsyi calculated using the different effective potentials 
is plotted, and the difference is shown to be very small. Thus, in this model the Landau 
gauge elementary field effective potential treatment and a PEP treatment give very similar 
numerical results at one loop level. 

We note that the similarities of the numerical results between the PEP and the Landau 
gauge fixed effective potential are not obvious a priori. In a gauge fixed formulation there 
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is a priori no reason to expect the results in one gauge to be numerically superior to those 
in another gauge. 

The similarities between the numerical results on the scale of new physics and therefore 
on the bound on the Higgs mass between the gauge invariant and gauge fixed formulation 
is due to the fact that the gauge couplings considered are relatively weak. For stronger 
coupling the gauge and scalar sectors will be more important in AA and the differences 
more noticeable. However, in that case a one-loop calculation will likely not be adequate 
and will thus be beyond the scope of our simple one-loop analysis. 

VII. CONCLUSIONS 

The usual formulation of bounds on the Higgs mass from the RG improved elementary 
field effective potential is afflicted with gauge dependence that renders it unsatisfactory. 
This gauge dependence is the result of the gauge dependence of the elementary field effec- 
tive potential itself, and it suggests that the elementary field effective potential is not the 
appropriate tool for the analysis. The error stemming from the gauge variance in calcula- 
tions in a specific gauge-fixing scheme cannot be inferred because the results can vary in a 
wide range by varying the gauge parameter. 

We have presented an alternative approach, using the recently-introduced PEP |T^, to 
calculate the Higgs mass bound in a simplified abelian model which displays the same vacuum 
stability problem as the Standard Model. We have explicitly provided a gauge invariant 
construction of the effective potential in terms of a gauge invariant order parameter that 
serves as a signal of symmetry breaking. This effective potential is unambiguously identified 
with the physical energy of a configuration and therefore provides reliable estimates from 
vacuum stability analysis. We constructed this gauge invariant effective potential to one 
loop order and improved it via RG. 

Compared to the analogous calculation in the Landau gauge the numerical results are 
similar when the gauge couplings are weak. While this suggests that Landau gauge cal- 
culations in the Standard Model may give numerically similar results to a gauge-invariant 
treatment, it does not provide justification for the principle of using the gauge-dependent 
effective potential to place bounds on gauge invariant quantities. 

Although not guaranteed a priori, our result on the numerical similarities between the 
scales and bound obtained from PEP and the Landau gauge-fixed effective potential suggest 
that in the case of weak gauge couplings, the Landau gauge effective potential provides 
a qualitatively reliable estimate. However, the PEP is necessary to establish the "error" 
estimate. In that respect, our results provide a tentative credibility to bounds based on 
Landau gauge effective potential calculations. 

An extension of the PEP formalism to nonabelian gauge theories would permit a gauge 
invariant calculation of the Higgs mass bound in the Standard Model using the PEP. Such 
a gauge invariant calculation is required to provide justified error bounds on the Higgs mass 
from vacuum stability considerations. We hope to provide such an extension in the near 
future. 



17 



ACKNOWLEDGMENTS 



D.Boyanovsky acknowledges support from the N.S.F. through grant No: PHY-9302534. 
We wish to thank Paresh Malde for useful correspondence and Tony Duncan for useful 
conversations. 



APPENDIX A: DIMENSIONALLY-REGULARIZED INTEGRALS 

Most of the integrals that arise in the calculation of the loop corrections to the gauge- 



invariant effective potential are familiar. The contributions of eqns. p.9| - |3.11| generate 
integrals of generic form 



(277 



m 



327r2 



1 ^ 3 



(Al) 



The more complicated term involving eq. p.l2| , however, requires additional effort to extract 
the 0{e^) piece. The angular integration gives 



■I (2^^ 



dk k 



d-2. 



2(27r)''-^ ""jo 



foo ^ 

/ du u yl + uyr + 
Jo 



u 



(A2) 
(A3) 



where fi^ 
written 



2^^772 



V[d] 



and r 



Ml 
'Ml 



Using eq.3. 197.1 and eq.9.111.2 of [27| this can be 



d-2 



(2vr: 



,d-l 





d 




'd-l 
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2-^1 




~2 




2 



1 d 3-d 



2' 2' 2 ' 

1 d-2 d + 1 



d-l 

+ r~T 



d-2 



d-l 



(A4) 



where d = A — e. The second term is finite for e ^ and so e = may be take immediately. 
It may be written as 



M^l 



477^4 



| — -y/r (1 + r) + (1 — r)^ arctanhi/rj . 



(A5) 



The first term contains a ^ pole arising from the F — | factor. A Laurent expansion about 
e = to 0{e^) is necessary, which requires the 0{e) terms of the hypergeometric function. 



This is obtained by writing 2F1 — — 2 + 



as a series expansion in r in terms of 



2' " ' 2' 2 ' ^ 

Pochhammer symbols, carrying out the expansion in e on each term, and resumming the 
resulting series. To C(e): 
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(4). 



9e 



e=0 



e=0 



(A6) 



Upon carrying through the differentiation the first term of eq.SE yields an infinite series 
that can be resummed. The second series terminates due to the factor (—2)^. The final 
result is: 



"2' ^2' 2 '\ 



:i-r) 



4r + (1 — r)^ In [1 — r] 



(AT) 



This is combined with the series expansions of the gamma functions and other prefactors, 
and terms to 0{e^) are collected. The result is: 



(27r 
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APPENDIX B: MS COUNTERTERMS AND /3 FUNCTIONS 
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2. (5 Functions 

/3a = (20A2 + 66* - 8!,* - Ue'X + i\y') (B7) 
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FIGURES 




0.5 

FIG. 1. log svi for the PEP as a function of e? and Xi^min for Vt = 0.5, k = v. 



DLog Svi - 




. 5 

FIG. 2. The difference between logsyi calculated from the PEP and logsyi calculated from 
the Landau gauge elementary field effective potential, as a function of ef and Xi^min, for Vt = 0.5, 

K = V. 
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